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Abstract 

We consider the matrix Aij, whose elements are independent ran- 
dom variables. We calculate the mean value of the number of the 
elements that we need to read to find miuj maxj A^j . 



1 The Min-max Algorithm 

In this paper we are going to analyze the process of finding minj maxj A^j 
[1,2,3], for which the reading of any single matrix element is particularly time 
consuming. However, this process-used in game theory and economy-may 
be shortened as to find the value of the Min-max, we do not have to read the 
whole matrix. Let us look at the following example: 

5 3 14 8 
2 4 9 1 
7 16 13 n 
12 6 10 15 

The maximal elements in each row are printed bold, and the elements which 
we do not have to read are underlined. 
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The algorithm of finding Min-max of matrix n x k reads 



inirmiax= infinity 
FOR i=l TO n 

mcLX= -infinity 

FOR j=l TO k 
t=A[i] [j] 

IF t > max THEN max=t 
IF t > minmax THEN BREAK 
IF max < minmax THEN minmax=max 

The essential instruction in this algorithm is: "IF t > minmax THEN BREAK" . 

The presented algorithm is a restricted version of algorithm alpha-beta 
pruning, that is used in logical computer games such as chess, Othello, etc. 

The search function of algorithm alpha-beta pruning reads [4,5] 

SEARCH (n, alpha) 

IF n = THEN RETURN - evaluation 
beta = - infinity 
FOR j=l TO m (list of moves) 
MOVE(j) 

t = - SEARCH ( n-1, - beta) 

UNDO MOVE 

IF t <= alpha THEN RETURN beta 
IF t < beta THEN beta=t 
RETURN beta 

2 Two equivalent models 

Generally, the time of performance of the algorithm depends on the distri- 
bution of matrix elements. 

In this paper we consider the probability model in which we assume that 
the matrix elements (^n, A„k) are independent uniformly distributed ran- 
dom variables on (0, 1). 

Our aim is to find the mean value of the number of matrix elements that 
are read while searching for the min-max. 
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We can also consider the discrete model in which the matrix elements 
{All, •••) Ank) are permutations of the set (1, 2, nk), each permutation be- 
ing of the same probability. 

For both models, the mean values are the same. 

Still another model is the discrete model in which the matrix elements 
belong to the finite set (i.e. the elements can be repeated). For this model, 
however, the obtained mean value is different from the ones achieved for the 
first two models. 



3 Recursion 

When the algorithm begins to operate, the parameter a takes the maximum 
value. After each row the parameter a can decrease. The further search 
depends on the value of a. Let M(a) stand for the mean value of the read el- 
ements, on condition that the algorithm begins to operate with the parameter 
a e (0,1). 

Now at the beginning of the matrix A^, let us add a new row {xi, X2, ...Xk) 
consisting of random independent variables uniformly distributed on (0, 1). 
Let us take an arbitrary number a G (0, 1). 

The probability that Xi,X2,Xj^i < a < xj is equal to (1 — ay^^. 

U xi,X2, ■ ■ ■ , Xj^i < a < Xj, the algorithm reads the elements xi, . . . , Xj, 
and, on average, M(a) elements from the remaining part of the matrix. 

The probability that {xi, . . . , x^} are in {t,t+A) is equal to (l + A)'=-i'=. 

If {xi, . . . , Xk} Q {t,t + A), the algorithm reads all k elements from the 
added row, and, on average, M{t) elements from the remaining part of the 
matrix. 

Therefore the mean number of the read matrix elements is equal to 



M{a) = (!-«)(! + M{a)) + (1 - a)(2 + M{a)) + ... + {l-a){k + M{a)) + 



pa (If 

+ (k + M(t))— dt^l + a + ... + a''-^ + M(a)+ 
Jo dt Jo 



dt 



dt. 



4 Calculations 

The obtained recursion makes it possible to calculate the mean number of 
the read matrix elements. 
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Let M„fc(a) stand for the mean number of the read matrix elements, on 
condition that the algorithm begins to operate with parameter a G (0, 1), n 
and k representing, respectively, the number of columns and number of rows 
in the matrix. 

Using the relation from the previous section we can write 

M„fe(a) = M„_i.fe(a) - T t^M„_i,fc(t) dt (1) 

Jo 

+1 + a + + ... + a''"^ 

This relation, together with the initial condition Mofc(a) = 0, allows us to 
find Mnk{0') for any n. 

Differentiating the equation ([T]) we obtain 

dM^ = (1 - ^.)«^zMM + (1 + 2a + 3a^ + ■ ■ ■ + (A: - l)a'~'). (2) 
da da 

We remember that 

dMr)u(a) , , , , 

^ 0, M„fc(0) = 0. (3) 



da 

Iterating n times the relation ([2]) we get 

= - a'y{l + 2a + 3a^ + --- + ik- l)a'-') (4) 

da 

= ^~^^r"'^" (l + 2a + 3a^ + --- + {k- l)a''-'). 
a'^ 

Hence 

Mnk{a) = r + 2t + + . . . + (A; - l)t''-^) dt = (5) 

JO t 



fc-l fc-ln-l , 

n 



1=0 1 = 1 j = 
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5 Results 



We obtain our main result, that is the formula for the mean number of the 
read matrix elements, by substituting a = 1 in the formula ([5]) 



k~l n—1 / \ 1 



Using the formula [6] 

«-)-"+e(,:0(-i>'-^- p) 



;i + a;) (l + a;)(2 + a;) (1 + x)(2 + - 1 + a;) ' 

we can rewrite our result in the following form: 

1=0 1=1 j=l r=l ^'^ ' 

For small values of n and m, we have 

M22(l)=3| M23(l) = 5^, M24(l) = 7^, 

1 31 241 Q 

M32(l) = 5-, M33(l)=7-, ^34(1) = 9^, 
M42(l) = 6§, M43(l)=9|?, M44(l) = 12-^^^^ 



35' 77' ^"^^ 45045 
Also, we would like to add that the following inequalities may be checked: 

Mnk{l) M„+i,fc(l) M„fc(l) M„,fc+i(l) 

n n + 1' /t fc + l' 

There is a good approximation for M„fc(l) 

n I n n \ 

^nfe(l) ~ 7T + ^ 1 + 



2 ylogn log n ) 

The error of the formula is less then 3 percent for > 28, > 17 , and less 
then 1 percent for n > 38, k > 20. 
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